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ABSTRACT 

A hydrostatic, primitive equation model with frontogenetical deformation 
forcing is used to simulate the passage of cold fronts over a two-dimensional ridge. 
The model includes a K-theory planetary boundary layer (PBL) parameterization 
with implicitly defined diffusion coefficients. Numerical simulations are performed 
for synoptic-scale ridges of varying widths both with and without frontal forcing. 
These results are compared to simulations that do not include a PBL 
parameterization, similar to previous inviscid studies by Williams et al. 

Relative to the inviscid results, the PBL simulations produced reduced 
frontolysis on the upwind slope and reduced frontogenesis on the lee slope, 
resulting in significantly smaller frontogenetic variations over the mountain. This 
is caused by convergence forcing in the well-mixed layer offsetting the overall 
frontolytical forcing on the upwind slope, and greatly reduced lee side convergence 
forcing due to the PBL. In contrast to the inviscid results, the final downstream 
front is weaker in the mountain simulations than in the flat-topography control case 
when PBL effects are included. In all PBL simulations, gravity wave generation 
is greatly reduced and no lee side hydraulic jumps are observed. In general, the 
inclusion of a PBL into the model results in more realistic wind and temperature 
fields compared to the inviscid model simulations. 
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I. INTRODUCTION 



Mountainous topography introduces dynamical complications to meteorological 
phenomena on all scales, from sub-mesoscale systems to the global circulation. Some 
effects such as diurnal mountain and valley winds are produced by the differential 
diabatic heating typically associated with mountainous terrain. Local wind regimes are 
also profoundly affected by sub-synoptic scale variations in topography. However, many 
important synoptic and mesoscale effects are related directly to the dynamics of 
atmospheric flow over synoptic-scale mountains. These effects include lee cyclogenesis, 
mountain wave generation, and downslope wind storms such as the Mistral in the 
southern French Alps. Equally important to forecasters is the effect of mountainous 
topography on atmospheric frontal intensity. Forecasters have long known that fronts tend 
to weaken as they move up the windward slope of a mountain range and often intensify 
on the lee slope. This effect has been documented in quantitative observational studies, 
for example in the analyses of high-resolution Alpine Experiment (ALPEX) data of 
Hartsough and Blumen (1990) and Radinovic (1986). 

Analytical solutions to the problem of fronts moving over mountain ridges are 
difficult because of the importance of ageostrophic circulations and the resulting 
nonlinearity of the equations of motion. Bannon (1983) derives frontal solutions for a 
quasi geostrophic front passing over a two-dimensional ridge. The linear Boussinesq 
equations he employs do not allow dynamic interaction with the mountain-forced 
circulations, but the results show windward slope weakening and lee slope strengthening 
of the cold front nonetheless. This is a result of the superposition of the mountain-forced 
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temperature field with that of the front. Bannon (1984) shows the effect of dynamic 
ageostrophic interactions in semigeostrophic solutions, but solutions are only obtained 
when the topography is allowed to vary in time. Blumen and Gross (1987) examine the 
effect of mountain dynamics by treating the front as passive scalar fields advected over 
a two-dimensional ridge by semigeostrophic steady-state fields. The scalar fields 
represent the temperature and tangential wind fields of a weak front. The solutions show 
the frontolytical effect of mountain-forced upwind divergence and the frontogenetical 
effect of lee side convergence.Other studies examine the dynamics of frontal interaction 
with topography by considering simple two-fluid models of the atmosphere (Davies 1984 
and Blumen 1992). 

While important realistic effects can be resolved by simplified analytical 
treatments, the full dynamics of frontal interaction with topography are analytically 
intractable. Furthermore, the simplified analytical solutions are only useful over the 
limited range of flow parameters for which the steady-state solutions are valid. Davies’ 
(1984) solutions, for example, are valid over a limited range of front-mountain aspect 
ratios. Pierrehumbert and Wyman (1985) characterize stratified flow over a ridge by the 
Rossby number (Ro) and a non-dimensional height parameter (h) of the mountain. They 
show that for small values of Ro and h, solutions based on the semigeostrophic 
approximations are valid. As h increases to near a critical height, however, the 
approximations break down and the flow exhibits vertically propagating gravity waves and 
a hydraulic jump on the lee side (Williams et al. 1992). Clearly, frontal solutions based 
on semigeostrophic flow are inaccurate for these flows. 
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Numerical simulation provides a means to investigate the dynamics of flows 
without analytical solutions. Zehnder and Bannon (1988) uses a semigeostrophic 
numerical model to study the effect of a front moving over a two-dimensional ridge 
toward a stationary downstream deformation field. The nonlinear mountain interactions 
are isolated and shown to produce a weakening frontal zone on the windward slope and 
intensification on the lee slope. This effect is attributed to temperature field deformation 
by mountain-forced convergence/divergence patterns. 

Williams et al. (1992) perform similar simulations using a two-dimensional 
primitive equation (PE) Boussinesq model. In these experiments a horizontal deformation 
field is applied to a thermal wind balanced disturbance superposed on a steady-state 
semigeostrophic mountain solution. The deformation field is allowed to move over the 
ridge with the mean background flow. The mountain-forced convergence/divergence 
fields produce frontogenetic variations similar to those found by Zehnder and Bannon. 
Both frontal and non-frontal simulations are performed over a range of mountain widths. 
The wider mountain non-frontal solutions vary little from the semigeostrophic initial state, 
but the steeper mountains produce vertically propagating gravity waves and lee side 
hydraulic jumps. This wave activity is amplified in the frontal simulations. Williams et 
al. also run the model without the deformation forcing to show the mountain effects on 
a passive tracer field similar to that considered by Blumen and Gross (1987). The passive 
temperature field gradient decreases on the windward side and increases on the lee side. 
In both the passive and forced frontal experiments, the mountain effect on frontal intensity 
is largely symmetrical (except for lee side gravity wave activity). In the passive 
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simulations, the final downstream intensity is nearly equal to the initial intensity and in 
the forced simulations, the final intensity in the mountain cases is similar to that of the 
flat topography case. 

A similar PE model is used by Gross (1994) to simulate a three-dimensional front 
passing over a finite, isolated synoptic-scale ridge. The front is produced by a developing 
nonlinear baroclinic wave moving over the ridge. In addition to demonstrating anti- 
cyclonic deformation of the frontal zone and horizontal flow diversion around the ridge, 
the study reveals three-dimensional effects relevant to cross-ridge flow. For low Froude 
number flows in which the flow is not blocked, the trailing synoptic anticyclone is 
impeded by the ridge, producing a cross-ridge pressure gradient. Down-gradient 
acceleration of the lee slope winds causes increased convergence and frontogenesis 
downwind of the ridge. Horizontal shear deformation is shown to have a small effect on 
frontal forcing over the ridge. Generally, however, the cross-mountain effects on the front 
agree well with the two-dimensional simulations of Williams et al. (1992). 

All of these studies consider inviscid flows only. While it is widely recognized 
that the planetary boundary layer (PEL) can have important effects on flow over 
topography, few studies address the problem quantitatively. Carruthers and Hunt (1990) 
review the effects of a uniform PEL capped by a strong inversion on the stratified 
overlying flow based on the linear analysis of Hunt et al. (1988) and the nonlinear 
analysis of Carruthers and Choularton (1982). Variations in the inversion height across 

the ridge are predicted, the nature of which depend on the Froude numbers at and above 

/ 

the inversion layer. The presence of a varying elevated inversion significantly affects the 
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mountain-forced flow accelerations within the PBL. The diffusive and momentum 



transfer effects of the PBL are not considered. 

Yang (1993) simulates two-dimensional stratified flow over a bell-shaped ridge 
using both hydrostatic and nonhydrostatic models with a turbulent kinetic energy (TKE) 
boundary layer parameterization. Large amplitude vertically propogating gravity waves 
are produced in both models and trapped lee waves are observed in the nonhydrostatic 
simulations. Abmpt hydraulic effects are not observed, however, except on extremely 
steep mesoscale ridges (slopes with height to length aspect ratios approaching 1.0). The 
study also shows that the hydrostatic and nonhydrostatic models produce similar flow 
stmctures when the the length scale of the ridge is much larger than a flow parameter 
defined as U/N, where U is the mean velocity scale and N is the Brunt-Vaisala frequency. 

The effect of the PBL has been treated in greater detail in flat topography frontal 
studies. Keyser and Anthes (1982) investigate PBL effects on frontogenesis using a two- 
dimensional hydrostatic PE model incorporating a mult-layer, first-order K-theory PBL 
parameterization. Detailed realistic frontal features not produced in inviscid simulations 
are evident in the PBL model results. The model is also run with a simple one-layer bulk 
drag PBL parameterization. The multi-layer K-theory PBL model is shown to produce 
more realistic results than either the bulk parameterization or the Ekman formulation used 
by Blumen (1980) in a similar model. Dunst and Rhodin (1990) obtain similar detail in 
their frontal simulations using a high-resolution first-order parameterization similar to the 
one used by Keyser and Anthes (1982), but with a different diffusion coefficient 
formulation. 
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This study seeks to improve the physical accuracy of previous numerical studies 
on frontal interaction with topography by including a realistic PBL parameterization. The 
basic model is the two-dimensional Boussinesq PE model of Williams et al. (1992, 
denoted W92) with the inclusion of a first order K-theory PBL parameterization following 
Keyser and Anthes (1982, denoted KA82). Numerous PBL closure schemes are evaluated 
by Holt and Raman (1988), who find that while the overall turbulence structure of the 
atmosphere is better modelled by more complex turbulent kinetic energy closure schemes, 
the mean sfructure of PBL is fairly insensitive to the closure scheme employed. The 
KA82 formulation is similar to the modified Djolov (1973) mixing length 
parameterization reviewed. In a separate study (publication pending), Peng and Williams 
compare the results of frontal simulations over flat topography using the KA82 
formulation to those using the Dunst and Rhodin (1990) approach. The two PBL 
parameterization s again produce very similar results. 

As in W92, the model is applied to infinite north-south ridges of varying widths. 
The simulations are run both with and without frontal deformation forcing in both inviscid 
and PBL modes. In most cases, model parameters are selected to correspond to those of 
W92 to facilitate comparison. The basic model and PBL formulations are developed in 
Chapter II. Model solutions are presented and discussed in Chapter HI, and the results 
are summarized in Chapter IV. 
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II. MODEL DEVELOPMENT 



A. BASIC EQUATIONS 

This study uses the W92 model modified by Peng and Williams to include a 
boundary layer formulation following KA82. Development of the inviscid equations in 
the following section follows W92. 

1. Inviscid Model Formulation 

The model is based on the adiabatic primitive equations with hydrostatic and 
Boussinesq approximations applied, which neglect the compressibility of the atmosphere. 
The domain is two-dimensional in x and z and is perpendicular to the semi-infinite ridge. 
Rotational effects are included in an /-plane approximation, and velocity and temperature 
fields are considered constant in y. This results in increasing errors as the equations are 
integrated in time, but W92 shows that neglecting along-ffont variations is a valid 
approximation in this frontogenesis study. The upper boundary is a rigid lid and the 
horizontal boundaries are cyclic. 

The equations of motion are modelled using a terrain-foUowing vertical coordinate 
C defined as 



z-z 




( 1 ) 



where z^ is the surface height, H is the height of the upper boundary and Z = H - z^ This 
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coordinate system places the lower boundary at C = 0 and the upper boundary at C = 1- 
The Boussinesq equations (Ogura and Philips 1962) are transformed with equation (1) to 
give 



dV • dV 2 

+ vvv + = -V(}) - -|-(C-1)0VZ -y^xv+F, 

ot dc 0 



( 2 ) 



60 • 60 

-.vve*c^=ft 



(3) 



3C 

V-(ZV) +7-^=0, 
dC 



(4) 



6(J) 




(5) 



where 



0 = 1U>jpf - 9 „, 



( 6 ) 






(7) 



In these equations, V is the horizontal gradient operator on C level surfaces. F and Q are 
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fourth-order diffusion terms for momentum and potential temperature, included to damp 
numerical noise in the inviscid model. The diffusion terms also balance the frontal 
forcing in steady-state to prevent sub-grid scale frontal collapse. They include terms of 
the form 




= -K 



d u 



4 ’ 
OX 



( 8 ) 



d^u 

F =-K , 

Vu 

dC 



(9) 



representing horizontal and vertical diffusion, with similar terms for v and 0. Q also 
represents the dry convective adjustment process. Coefficients and Ky are given 
simple vertical profiles in the inviscid model. Other symbols have the usual 
meteorological meaning. Vertical boundary conditions are 

C = 0 applied at C = 0, 1. 



Horizontal boundary conditions are periodic. 

The frontogenesis simulations are forced by the vertically uniform horizontal wind 
deformation field 



= (D/j.Ocosh(t^)cos[p(x - Ut) + Xj], 
= (Z)/^J$sinh(^^)sm[^<J: - Ut) + 



( 11 ) 
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where D, ju = 2nAV, and x, are constants and U is the mean background flow. This 
deformation field is cyclic in x and moves with the mean flow. It is also nondivergent 
and irrotational. 

The dependent variables can be expanded into mean, deformation, and perturbation 
components as follows: 

\ = {a[U^UJ.x- Ut,y)] + u(xX,t)} i 

^{aVJix-Ut,y)^vix,C^)}h 

C = C(:c,C,0. 

0=0(x,CA 

4) = - Ut,y) -fUy + (p(x,C,f), 

where 

a = (Z)/Z. (13) 

The brackets (( )) indicate a horizontal average over the domain and a is factor accounting 
for the effect of atmospheric thickness variations on the mean and deformation flows. 
The deformation pressure function maintains the geostrophic balance of the 
deformation flow. 

Equations (2) - (5) can be put in flux form by multiplying equations (2) and (3) 
by Z and substituting from equation (4), using the vector identity 

V- (/4 V) = V + (Vy4)- V, (14) 
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where A is a scalar. Substituting the variable forms of equations (12) and evaluating at 
y = 0 yields 



+ 4-[Z(at/+«)^ + 4r[Z(aI/+«)a + [a f/ fa !/+«)] 
dt dx dC dx ^ 

3(p g(C-l)0 dZ 

-ZaU +Z—- ^\ Z—-PZ-ZF =0, 

dx dx 0 dx “ 



(15) 



dv d 0 1 0v 

^ ^ + Za U — 

dt dx d(, *^dx 

d 

- V— (Za U^-ZfU ^J{aU*u)Z- ZF^ = 0, 



(16) 



00 a a • 00 

Z — + — [Z(a U+u)d] + — (Z0C) +ZaU — ~ZQ =0, 
dt dx aC ‘‘dx 



(17) 



^[Z(af/+«)]+z|^=0, 

0x ac 



(18) 





(19) 



The pressure perturbation (p must be eliminated to close this system of equations. 



Integrating equation (19) over C yields 
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( 20 ) 






where 



C 




( 21 ) 



and (Po is the surface value. Equation (18) is then integrated from C = 0 to C = 1 to give 



where the overbar indicates vertical integration over the entire domain. To obtain a 
closed form we also assume that the mass flux does not vary in time so that 



This assumption neglects the effect of mountain wave drag and the loss of mean flow 
momentum to the boundary layer. These minor losses are assumed to be offset by large- 
scale forcing and are not explicitly included in the model. 

The pressure function form of equation (20) is introduced into equation (15) to 

give 



dx 



d — 
—(uZ) = 0, 



( 22 ) 



(ttZ) = 



(23) 




^<pQ ^ g(C-i)0 ^az 



dx dx 0 dx 
0 



Z — ^fvZ^ZF , 
Fit “ 



(24) 



dx 
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where L contains the momentum flux terms in equation (15). Finally, we apply the 
constant mass flux assumption of equation (23) to the vertical average of equation (24) 
and subtract the result from equation (24), giving 
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(25) 



e dZ _ 

-|-[(C - 1)0 - (C - 1)0]Z— - V)Z + Z(F -F ). 

O ox u u 



0 



Equation (25) provides a prognostic equation for u in terms of V and 0. The model is 
based on the numerical integration of equations (25), (16) and (17) in time. 

2. Boundary Layer Formulation 

The boundary layer parameterization employed follows KA82, which is based on 
the high-resolution nocturnal PBL parameterization of Blackadar (1978). Turbulent fluxes 
are represented in K theory by 




(26) 




(27) 
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( 28 ) 



-7T/ 50 

w'Q'=-K—. 

dz 



The diffusion coefficients and are calculated implicitly from model shear and 
stability fields to minimize a priori assumptions about the boundary layer structure. 
Blackadar (1978) derives a closed system of equations for calculating and Kq. based 
on the second-order closure scheme of MeUor and Yamada (1974). He adapts the Mellor 
and Yamada “level-2" approximation, which neglects advective and diffusive terms in the 
second moment turbulence equations, such that the equations obey Monin-Obukov 
similarity. The equations are solved in terms of the mixing length /, vertical shear s, and 
the Richardson number Ri, given by 



Ri = ||^T 
eaz 



S = 



dV 

dz 



( 29 ) 



In this formulation, / is not explicitly calculated and must be specified based on empirical 
evidence of the mixing length above the surface layer. Over the range of positive and 
slightly negative Ri, the resulting distribution is approximated by 



l.l(/^5) 



K 

mz 




= < 



Ri -Ril 

C 




0 , 



Ri<Ri 

I 

Ri^Ri 



( 30 ) 
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where the critical Richardson number Rij = 0.25 (Blackadar, 1978). Note that the vertical 
shear dependence is included in Ri. Following KA82, the coefficients are broken down 
into a small, constant diffusive part and a variable part modelled after equation (30): 

{ Ri -Ri'l 



K 

mz 



0z 



K 



zO 



I 2f2 

I S 

0 




Ri<Ri 

c 



( 31 ) 




RikRi 

C 



where kQ is the von Karman constant. 

This formulation defines the vertical fluxes of heat and momentum above the 
surface layer, which is assumed to be contained within the first model layer. The surface 
layer fluxes are modelled after Monin-Obukov similarity theory, which gives 



/ 

u w 



/ 



/ 

V w 



/ 



2 

u u 



/ 2 2 . 1 / 2 ’ 
(M +Vl 

2 

U V 






( 32 ) 



w^0^ = O, 



where 



ln(z/z)-t|r 

u ffi 



( 33 ) 



( 34 ) 
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Here Zo is the surface roughness length (the height at which velocity decreases to zero) 
and i|j„ is the surface layer stability correction term. Diabatic surface heating is not 
included in the model, giving the lower boundary condition of equation (33) and requiring 
a neutral surface layer so that i)j„, = 0. Other constants are listed in Table 1. In equations 
(32) through (34) u and v are the velocities at the first model level and z is the 
corresponding height. 

Given the turbulent fluxes computed from equations (26) through (28) and (32) 
through (34), the basic momentum and heat equations are modified by adding the terms 



diu'w') 



1 - 



acv'w') . 

ac 



and 



3(w^0^) 



(35) 



to the right hand sides of equations (2) and (3), respectively. The fourth-order diffusion 
coefficients in equations (8) and (9) are also modified as described in the next section. 
B. COMPUTATIONAL ASPECTS 

The model equations are solved numerically using finite differencing on an 
Arakawa B grid (Arakawa and Lamb, 1977). Temperature, pressure, and vertical velocity 
points are offset horizontally from horizontal momentum points, and vertical velocity 
points are also offset vertically. Finite differences are centered in space and time, with 
a Matsuno scheme (Euler backward) time step inserted every fourth iteration to control 
solution separation. Convective adjustment is applied every three time steps by vertically 
averaging the potential temperature field in areas of static instability. The model domain 
extends 3600 km in the east- west (x) direction and 12 km vertically. Horizontal grid 
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Table 1. Specification of constant parameters. 



Parameter 


Numerical value 


Parameter 


Numerical value j 


W 


3600 km 


Rio 


1 


H 


12 km 


h 


0.4 I 


Ax 


40 km 


Ko 


1.4 X 10" m'' s ' 


AC 


1/50 


Kho 


10^ m^ s' 


At 


90 s 




1 m' s-' 


h 


2 km 


Zo 


0.4 m 


U 


10 m s' 


1 


100 m 


D 


10-^ s ' 


f 


10-^ s ' 


BQJdz 


4 Kkm-' 


g 


9.8 m s" 


a 


12 K 




1003 m^ s^K' 


p 


1.75 X 10-^ km ' 


00 


300 K 


Xi 


- TC 


Po 


1000 mb 


Ct 


0.75 





spacing is 40 km and there are 50 vertical levels uniformly spaced in z, corresponding to 
240 m spacing over flat topography. The time increment is 90 s. Other numerical 
constants are listed in Table 1. The surface topography is defined by 



. 2 
z =hcos 

5 



nix-W/2) 



\x-W/2\ 



eW 




(36) 



The mountain height h is held constant while the width e (the fraction of the horizontal 
domain occupied by the mountain) is varied between 0.6 and 0.2. 
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Fourth-order diffusion constants in the inviscid model are given the value Kq at 
the surface and increase linearly in the vertical to a factor of 10 at the upper boundary. 
The linear increase is included to maximize the damping of small-scale numerical noise 
throughout the domain without including excessive surface diffusion. A hyperbolic 
sponge layer is also included near the top of the domain. The output fields show little 
damping of resolved waves and the sponge layer effectively eliminates reflection off the 
upper boundary. In the PEL model, the fourth-order diffusion coefficients are defined 
following KA82 as 






r 

1 2 * 1 


'3u' 


2 


^ dv^ 


2 


1/2 


K + 

HO 2 ® 




+ 






► 



(37) 



Vertical and horizontal dilfusion coefficients for u, v and 0 are given equal values within 
each model. 

Boundary layer diffusion is computed using equations (31) through (35), with the 
constant values given in Table 1. Finite differencing tends to give systematically greater 
values of Ri than the mean Ri within the layer (Blackadar, 1978), so Ric = 1.0 vice the 
theoretical value of 0.25 to account for the finite grid resolution (following KA82). 

Initial fields are the same as those used in W92. The non-frontal mountain initial 
flow is defined by semigeostrophic solutions (following Merkine, 1975) for constant U 
and stratification away from the ridge and periodic boundary conditions. In frontal 
simulations, a perturbation is added to these basic fields, defined by 
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( 38 ) 



Here dQjdz is the basic field stratification and a is the perturbation amplitude. Note that 
the perturbation is limited to the lower part of the domain (below Cr) to prevent 
frontogenesis at the upper boundary. It is also offset horizontally such that the 
maximum perturbation coincides with the axis of dilatation of the deformation field 
(equation (11)). Initial u, v and w fields are derived using the thermal wind, 
quasigeostrophic circulation, and continuity equations. Details are given in W92. 
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III. NUMERICAL SOLUTIONS 



The results of the numerical simulations are presented and discussed in the 
following sections. In each case, the inviscid data are similar to those presented in W92, 
which provides a complete discussion of these results. The emphasis of this study is the 
comparison of the inviscid case results to those of the PEL simulations, which are 
integrated from the same initial conditions. 

A. FRONTAL SOLUTIONS WITH NO TOPOGRAPHY 

The model is run first with no topography to serve as a control case for the 
mountain simulations. The initial fields are shown in Figure 1 as functions of x and z. 
The cross-frontal velocity field in Figure lb (and all subsequent figures) includes the 
mean background flow U but not the constant deformation field component. At time t 
= 0, the axis of dilatation of the deformation field is coincident with the maximum 
surface temperature gradient at x = 900 km. The stretching deformation, which moves 
with the background flow, will therefore intensify this cold frontal zone. 

Figure 2 contains the frontal solutions after 48 hours of integration. Only the 
lower part of the atmosphere is shown to highlight near-surface effects. The front has 
intensified and all fields have been advected downstream with the background current. 
These fields cleariy show the maximum thermal gradient and associated vorticity maxima 
at the surface, decreasing rapidly with height. The corresponding PEL solutions are 
shown in Figure 3. Strong vertical mixing is evident in the near-vertical potential 



temperature contours throughout the 1.0 to 1.2 km well-mixed layer. The frontal zone 
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(d) potential temperature, deg C. 
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V at t = 0 b. U at t 



V at I = 48 b. U at t = 48 
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Figure 2. Inviscid frontal solutions at t = 48 h. 
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Figure 3. PBL frontal solutions at t = 48 h. 
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V at t = 48 b. U at I = 48 




is weaker than that shown in Figure 2, and the thermal gradient is constant throughout the 
PEL. Here the cross-front and along-front wind maxima are elevated jets near the top of 
the PEL as winds are frictionally forced to zero at the surface. 

In a separate study (discussed in Chapter I), Peng and Williams obtain frontal 
solutions similar to those shown in Figure 3. Employing a variable horizontal coordinate, 
they present a more detailed depiction of the frontal structure. A major finding is that 
the simulated fronts are always weaker when PEL effects are included, due primarily to 
vertical mixing. The turbulent mixing of the weaker-gradient upper PEL air reduces the 
potential temperature gradient at the surface. This effect can be seen in Figures 2 and 3. 
The PEL also produces diffusion in addition to the fourth-order diffusion, resulting in 
further smoothing of the fields and a reduction of the thermal gradient. 

To quantify the strength of the frontal zone, the parameter d is defined as 

. |A9| 

d= , 

00 (39) 

dx EDiaX 

where A0 is the maximum horizontal potential temperature variation on the lowest 
numerical level. The temporal evolution of the value for the flat topography 
simulations is shown in Figure 4. After an initial period of reduced frontogenesis, the 
PEL front undergoes frontogenesis at a rate nearly equal to that of the inviscid front 
(indicated by the slopes of the curves in Figure 4) until about 36 h, when it reaches 
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maximum intensity. At this point the frontogenetical forcing is matched by the diffusive 
forcing and the front is in steady-state. The curves in Figure 4 will be used as a 
benchmark for the mountain simulations. 

The front in the PEL case also translates slower than the frictionless front, which 
is over 200 km further downstream at t = 48 h. The difference between the frontal 
positions (defined by the location of maximum potential temperature gradient) is shown 
as a function of time in Figure 5a. The increase in the distance that the PEL front lags 
behind is correlated to the difference in maximum cross-frontal velocities (Figure 5b). 
The velocities shown are vertical averages over the four lowest numerical levels 
(equivalent to 960 m) and are representative of the average PEL cross-frontal wind at the 
frontal surface. The fronts are separated by no more than 50 km through t = 18 h, when 
the velocities are similar. After t = 18 h, the velocities diverge and the frontal separation 
steadily increases. Thus it appears that the reduction of the cross-frontal wind velocity 
caused by the PEL results in a slower translation speed and a different vertical tilt. 
Eecause the deformation field moves with the background flow, however, it is expected 
that both fronts will eventually move with the same speed. Clearly, vertical mixing in 
the PEL plays a role in the causing the PEL surface front to lag behind the inviscid front. 
B. NON-FRONTAL SOLUTIONS WITH TOPOGRAPHY 

Eefore investigating topographic effects on the fronts, the effects of the PEL on 
the basic flow over the ridge must be examined. In these experiments, the 
semigeostrophic initial fields are integrated in time with no temperature perturbation or 
deformation forcing applied. The initial fields for the e = 0.6 ridge are shown in 
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Frontal Strength 




Figure 4. Frontal strength parameter for flat topography; solid line: inviscid case; 
dashed line: PEL case. 
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a. Difference in Frontal Position 




b. Maximum Cross-Frontal Velocity 




Figure 5. Comparison of frontal separation with cross-frontal velocity in flat 
topography cases: (a) inviscid frontal position minus PBL frontal position, km; 

(b) maximum cross-frontal velocity, m/s. Solid line: inviscid case; dashed line: PBL 
case. 
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Figure 6. The w contours in Figure 6c are of true vertical velocity in height coordinates. 
The fields are symmetrical about the ridge axis, as expected. The t = 48 h fields (Figure 
7) show only a slight departure from the semigeostrophic solutions as weak wave activity 
is generated and the cross-ridge velocity maximum shifts slightly toward the lee slope. 
The PEL solutions for t = 48 h (Figure 8) are more asymmetrical in the lower layers in 
this respect as the elevated cross-frontal jet is displaced significantly over the lee slope. 
No gravity waves are evident in the figure, however. The 6 = 0.4 ridge produces similar 
t = 48 h solutions (Figures 9 and 10), but with stronger mountain-forced wind fields and 
larger amplitude gravity waves, which begin to appear in the PEL case (Figure lOd). 

The e = 0.2 mountain inviscid solution exhibits a dramatic departure from the 
semigeostrophic solution. The fields at t = 48 h are shown in Figure 11. As the 
mountain slope becomes steeper, the cross-mountain flow is accelerated more strongly and 
advected over the lee slope. In the inviscid simulations, the maximum winds are at the 
surface and an evident hydraulic jump is generated when the high speed winds reach the 
base of the slope. This causes the well-mixed area at the base of the lee slope and the 
flow reversal above the lee slope u and v maxima. Large amplitude, upward propagating 
gravity waves also result from the hydraulic jump. In contrast, the PEL solution at t = 
48 h (Figure 12) shows no such abrupt effects. Lee slope winds are enhanced, creating 
gravity waves, but the elevated jet is “insulated” from the topography by the PEL, so the 
topographic slope changes do not force the abrupt changes in momentum exhibited in the 
inviscid simulation. 
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Figure 6. Initial fields for non-frontal mountain (€ = 0.6) case. 
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V at t = 0 b. U at t 
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Figure 7. Inviscid non-frontal mountain (e = 0.6) solutions at t = 48 h. 
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Figure 8. PBL non-frontal mountain (e = 0.6) solutions at t = 48 h. 
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V at t = 48 b. U at t = 48 




V at t = 48 b. U at t = 48 
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Figure 9. Same as Figure 7 except for e = 0.4 mountain. 
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Figure 10. Same as Figure 8 except for e = 0.4 mountain. 
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V at t = 48 b. U at t = 48 




V at t = 48 b. U at I = 48 
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Figure 11. Same as Figure 7 except for € = 0.2 mountain. 
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Figure 12. Same as Figure 8 except for e = 0.2 mountain. 
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V at t = 48 b. U at t = 48 





In order to investigate the dynamical effects of the mountain-forced circulations 
on frontogenesis, the divergence fields are examined. The divergence fields of the e = 
0.6 mountain solutions at t = 24 h are shown in Figure 13. The mountain circulations are 
weU-developed by 24 h, and this roughly corresponds to the time that the fronts cross the 
ridge in the frontal simulations. In the inviscid case (Figure 13a), the flow is divergent 
on the upwind slope as the air accelerates over the ridge. Strong convergence is shown 
on the lee slope, caused by the rapid decrease in velocity downstream of the cross-ridge 
jet velocity maximum. The horizontal divergence field is complicated in the PBL case 
(Figure 13b) because the wind decreases to zero at the surface, producing significant 
convergence in the PBL on the upwind slope overridden by a divergent area as the winds 
above the PBL accelerate over the ridge. This pattern is consistent with the vertical 
velocity pattern, which increases from zero at the surface to a maximum just above the 
PBL. Conversely, the lee slope is characterized by divergence within the PBL overridden 
by an area of convergence. There is less convergence on the lee slope in the PBL case 
because the winds do not decrease so rapidly near the base. The weaker elevated jet does 
not enhance the lee side convergence surface as much as in the frictionless case. 

The boundary layer also forces a horizontal temperamre gradient within the PBL 
due to vertical mixing on the mountain slope. Vertical mixing reaches a higher 0-level 
near the top of the ridge than near the base because the stability increases over the ridge. 
The PBL is therefore mixed with air at a higher potential temperature near the top of the 
ridge. This produces a positive temperature gradient on the upwind slope and a negative 
temperature gradient on the lee slope, as shown in Figure 14c at t = 12 h. The thermal 
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a. Divergence (Inviscid) at t = 24 




b. Divergence (PBL) at t = 24 




Figure 13. Horizontal divergence at t = 24 h for non-frontal mountain (e = 0.6) 
cases: (a) inviscid case, 10'^ s'^; (b) PBL case, 10'^ s'^. 
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Figure 14. Potential temperature fields for non-frontal mountain (e = 0.6) cases: (a) inviscid case at t = 12 h; (b) inviscid 
case at t = 24 hr; (c) PBL case at t = 12 h; (d) PBL case at t = 24 hr. 



pattern advects downstream somewhat, but is nearly at steady-state by t = 24 h (Figure 
14d). The potential temperature fields for the same times in the inviscid case show no 
significant horizontal temperature gradient along the mountain surface (Figure 14a-b). 
The effect is unique to the mountain PBL cases because mixing to differing 0-levels is 
required to produce a horizontal temperature gradient. The strength of the gradient 
produced is also highly dependent on the vertical thermal structure, as will be 
demonstrated in the frontal simulations. 

C. FRONTAL SOLUTIONS WITH TOPOGRAPHY 

1. Intensity Variations 

In these simulations, the frontal forcing is applied to the model over the mountain 
domain. The initial conditions are a superposition of the mountain semigeostrophic 
solutions (Figure 6) and the frontal initial conditions shown in Figure 1. In all cases the 
axis of the frontal forcing begins at x = 900 km. Figures 15 though 17 contain the time 
evolution of the potential temperature field at 6 h intervals for both the inviscid and PBL 
frontal simulations over the 6 = 0.6 mountain. The initial temperature perturbation 
(Figure 15a and 15d) has a positive gradient at the surface on the western half of the 
domain and a negative gradient on the eastern half. By t = 6 h, a well-mixed layer is 
already developed in the PBL case (Figure 15e). The vertical mixing has produced a 
significant strengthening of the surface potential temperature gradient compared to the 
frictionless potential temperature field (Figure 15b). Here the mixing has a greater effect 
on the thermal gradient than in the non-frontal case because the initial temperature 
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a. Theta (Inviscid) at t = 0 



d. Theta (PBL) at t = 0 




b. Theta (Inviscid) at t = 6 






c. Theta (Inviscid) at t = 12 f. Theta (PBL) at t = 12 





Figure 15. Potential temperature fields for frontal mountain (e = 0.6) cases: 
(a)-(c) inviscid case for t = 0 to t = 12 h; (d)-(f) PBL case for t = 0 to t = 12 h. 
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a. Theta (Inviscid) at t = 18 



d. Theta(PBL)att= 18 




b. Theta (Inviscid) at t = 24 




c. Theta (Inviscid) at t = 30 





e. Theta (PBL) at t = 24 




f. Theta (PBL) at t = 30 




Figure 16. Same as Figure 15 except for time t = 18 through t = 30 h. 
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Figure 17. 



d. Theta (PBL) at t = 36 




e. Theta (PBL) at t = 42 




f. Theta (PBL) at t = 48 




Same as Figure 15 except for time t = 36 through t = 48 h. 
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perturbation increases the differential 0-level mixing effect described in the previous 
section. The trend continues through t = 12 h as the well-mixed layer continues to 
develop and deepen. 

At about t = 18 hr, the inviscid front reaches the crest of the ridge and moves 
down the lee slope in the next 18 h, undergoing rapid frontogenesis (Figures 16a-c and 
17a). The PEL front lags behind by about 12 h, and increases in intensity at a much 
lower rate as it reaches the lee slope (Figure 16d-f) and then gradually weakens on the 
lower half of the lee slope (Figure 17d-f). Conversely, the inviscid front continues to 
strengthen over the entire lee slope and weakens slightly as it reaches the downwind plane 
at t = 48 h (Figure 17c). By the end of the simulation at t = 48 h, the PEL front is much 
less intense than the inviscid front. 

The frontal strength variations are summarized for the e = 0.6 mountain in Figure 
18, which contains the temporal evolution of the values for both frictionless and PEL 
simulations. The flat topography frontal solutions are included in the figure for reference. 
The corresponding six-hourly frontal positions are shown in Figure 19. The values 
show that in the frictionless case, the front weakens significantly on the upper half of the 
upwind slope and then undergoes strong frontogenesis on the lee slope. This is the same 
effect observed in W92, attributed to the mountain-forced divergence field. The 
weakening near t = 48 h is also observed in W92 and is shown to be caused by a small 
area of divergence at the base of the lee slope. The PEL rf-value curve confirms that the 
PEL front is stronger on the upwind slope and weaker on the lee slope compared to the 
frictionless case. 
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Frontal Strength 




Figure 18. Frontal strength parameter for 6 = 0.6 mountain; solid line: inviscid flat 
topography case; dashed line: PEL flat topography case; dot-dashed line: inviscid 
mountain case; dotted line: PEL mountain case. 
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a. Frontal Position, Inviscid 




b. Frontal Position, PBL 
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Figure 19. Six-hourly frontal positions corresponding to asterisks in Figure 18: (a) 
inviscid case; (b) PBL case. 



Comparison of the d-value curves to those of the flat topography simulations 
reveals several interesting features. In the flat topography cases, the PBL front is weaker 
than the inviscid front at all times. Note also that the inviscid front is weaker in the 
mountain case than in the flat topography case even on the lower half of the upwind slope 
(before the period of frontolysis). The mountain PBL front, however, is stronger than 
both the inviscid mountain front and the flat topography PBL front for the first 12 h of 
the simulation. Thus the combination of the mountainous topography and the PBL has 
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a frontogenetical effect on the upwind slope. The overall effect opposes this upwind 
slope effect, however. The final strength of the fronts in the frictionless simulations are 
similar whether or not topography is included (confirmed in W92 by extending the 
integration to t = 60 h). In the PEL simulations, however, the front is significantly 
weaker in the mountain case at t = 48 h. In contrast to the inviscid results, when PEL 
effects are included, the mountain appears to have a net frontolytical effect. 

The steeper mountain (e = 0.4) if- value plot (Figures 20 and 21) confirms the 
above analysis. Here the inviscid front is stronger on the upwind slope, but this is due 
to frontogenesis prior to and right at the base of the ridge. This is caused by a 
convergence area near the base of the upwind slope (as found in W92). The PEL front 
is less affected because the PEL effect forces the frontal zone immediately on to the 
mountain slope (Figure 21). Inspection of the slopes of the curves reveals that the PEL 
front undergoes greater intensification and less weakening than the inviscid front once 
they are over the upwind slope. The e = 0.2 mountain simulation produces similar trends 
(Figures 22 and 23). Detailed analysis of the mountain-forced intensity variations is not 
possible with six-hourly time resolution in this case because the fronts move over the bulk 
of the ridge in 6 h or less. 
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Frontal Strength 




Figure 20. Same as Figure 18 except for e = 0.4 mountain. 
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Figure 21. Same as Figure 19 except for e = 0.4 mountain. 
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Frontal Strength 




Figure 22. Same as Figure 18 except for e = 0.2 mountain. 
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Figure 23. Same as Figure 19 except for e = 0.2 mountain. 
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2. Frontogenetical Forcing 

The dynamical mechanisms responsible for the observed variations in frontal 
intensity can be determined by examining the terms of the frontogenetical forcing 
function, which can be written as 
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The convergence term (equation (42)) does not include the confluence of the deformation 
field, which is given by the deformation term (equation (44)). The form of equation (45) 
is from KA82. The terms represented include forcing due to PEL, numerical fourth-order 
diffusion, and convective adjustment diffusive processes. These terms are not calculated, 
so a complete summation of forcing terms cannot be computed to explain frontal intensity 
variations. Rather, the terms and are used to determine the differences in 

forcing between the frictionless and PEL simulations. The figures presented in this 
section do not include the advective term (equation (41)) and should be considered 
Lagrangian forcing fields (acting along parcel trajectories). All terms are calculated in 
z coordinates and the forcing fields in the figures have been filtered to remove grid-scale 
noise introduced in the numerical derivatives of the interpolated fields. 

The frontogenetical forcing fields for the e = 0.6 ridge at t = 6 h are shown in 
Figure 24. The PEL simulation produces relatively large convergence forcing within the 
boundary layer on the upwind slope (Figure 24d). This is caused by the convergence 
within the PEL and compounded by the stronger thermal gradient induced by vertical 
mixing (note the 60/3x factor in equation (42). The convergence term in the inviscid 
simulation shows weak frontolytical forcing (Figure 24a) due to divergence. The 
compounding effect of the PEL temperature gradient is also observed in the deformation 
forcing field (Figure 24c and 24f). Here the confluence in the PEL case is equal to that 
in the inviscid case, so the increased forcing in the PEL case is caused solely by the 
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Figure 24. Frontogenetical forcing fields at t = 6 h for e = 0.6 mountain: 

(a)-(c) inviscid case; (d)-(f) PBL case (all in 10’‘° deg C m'‘ s'*). Dotted lines are 
potential temperature contours; an asterisk indicates the location of the surface front. 
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greater temperature gradient within the PBL. The tilting term is frontolytical in both 
cases (Figure 24b and 24e) on the upwind slope and is larger in the PBL simulation, 
partially offsetting the increased convergence forcing. 

Figure 25 contains the forcing fields at t = 18 h, when the fronts are near the crest 
of the ridge. The inviscid front is weakening because of the negative convergence and 
tilting forcing upwind of the ridge crest. At the top of the ridge, the tilting term is 
strongly frontogenetic and dominates the convergence forcing (Figure 25a-b). The PBL 
smoothes the transition from upward to downward vertical motion at the crest (reduces 
dwidx in equation (43)) and produces a weaker tilting effect. Forcing is nearly balanced at 
the PBL front as indicated by the neutral trend at t = 18 h in Figure 18. Convergence 
forcing is still positive but weakening near the crest of the ridge, and tilting forcing is 
much weaker than in the inviscid case (Figure 25d-e). The tilting term is positive at all 
times over the top of the ridge and in fact does not vary much in magnitude as the front 
moves over the ridge. This suggests that the tilting effect is produced largely by the basic 
flow over the ridge. The forcing fields for the basic flow at t = 18 h are shown in Figure 
26. Comparison of the tilting term field at the top of the ridge shows little difference in 
the magnitude of the forcing produced by the basic field alone and that of the frontal 
case. 

As the fronts move on to the lee slope (Figure 27, at t = 24 h), the inviscid front 
enters the strong convergence area producing strong positive convergence forcing. In 
agreement with Gross (1994), tilting forcing becomes significantly frontolytical as the 
front passes through the lee side vertical velocity maximum, but in this case the positive 
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Figure 25. Same as Figure 24 except for t = 18 h. 
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Figure 26. Same as Figure 25 except for non-frontal cases. 
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Figure 27. Same as Figure 24 except for t = 24 h. 
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convergence forcing is nearly twice as large in magnitude (Figure 27a-b). In the PEL 
case, convergence forcing is weak and negative within the PEL (due to lee slope PEL 
divergence) and positive in the convergence area near the top of the well-mixed layer 
(Figure 27d). The well-mixed layer extends through both the positive and negative 
forcing areas, so the effect of the convergence forcing field is unclear. The tilting effect, 
for example, is nearly neutral within the boundary layer at all times, but affects the 
surface front by forcing the potential temperature gradient at the top of the PEL. It is 
suspected that the PEL forcing term (equation (45)) counteracts the convergence forcing 
to some degree within the PEL. Overall, the convergence forcing in the PEL simulation 
is dramatically weaker on the lee slope than in the inviscid case. 

The forcing fields at t = 42 h are contained in Figure 28. The inviscid fields are 
complicated because of the gravity waves generated when the front reaches the base of 
the lee slope, but the forcing is nearly balanced at the front because the intensity is not 
changing in time at t = 42 h (Figure 18). The PEL front is undergoing a gradual 
weakening trend from t = 36 h to t = 48 h. Strong frontolytic tilting forcing is present 
over the lower half of the lee slope (Figure 28d). The magnitude of the convergence 
forcing (again negative near the surface and positive near the PEL top) is considerably 
less than that of the tilting term. Deformation forcing is positive but a factor of ten 
smaller in magnitude. Thus the strong negative tilting effect is not balanced by positive 
convergence forcing as it is in the inviscid case and the PEL front weakens. Diffusive 
forcing is also probably a significant frontolytical factor in this region because the other 
forcing is relatively weak. 
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Figure 28. Same as Figure 24 except for t = 42 h. 
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rv. SUMMARY OF RESULTS 



The numerical simulations contained in this study demonstrate that the PBL plays 
an important role in frontal dynamics over mountainous topography. The effects in the 
flat topography case are consistent with those found in similar studies. In the mountain 
simulations, the elevation of the cross-ridge jet to the top of the PBL greatly reduces 
gravity wave generation and eliminates hydraulic Jump effects in the mountain profiles 
considered. These results are more realistic than those of the W92 inviscid study, which 
generates more wave activity than is realistic for smooth, synoptic-scale ridges. The 
weaker elevated jet also results in much weaker lee slope convergence, which is the 
dominant forcing mechanism in the inviscid study. Turbulent mixing on the mountain 
slopes is shown to indirectly affect frontal dynamics. The PBL reaches higher 9-levels 
on the ridge than on the adjacent planes, producing a horizontal temperature gradient in 
the well-mixed layer. This gradient has a large effect on frontal forcing because the 
convergence terms are proportional to the horizontal temperature gradient. 

In summary, the front still shows ffontolysis on the upwind slope and increased 
frontogenesis on the lee slope when PBL effects are included, but the magnitude of the 
variation is less than in the inviscid simulations. On the upwind slope, frictional 
convergence in the PBL due to the mountain slope combines with the mixing-enhanced 
horizontal temperature gradient. This produces a stronger front on the upwind slope and 
partially counteracts the weakening caused by the divergence associated with the 
accelerating cross-mountain flow above the PBL. At the crest of the ridge. 
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frontogenetical tilting is the dominant forcing mechanism and is weakened by the addition 
of the PBL. On the lee slope, convergence forcing is quite small or even negative in the 
PBL simulations because of reduced convergence associated with the cross-ridge jet and 
frictional divergence in the PBL near the surface. The result is that the PBL front 
intensifies much more slowly than in the inviscid case and begins to weaken over the 
lower half of the lee slope. In contrast to the inviscid simulations, the net effect of the 
front passing over the mountain is ffontolytical when PBL effects are included. Narrower 
mountain profiles produce slightly different results because the fronts begin forming 
further upstream relative to the mountain, but the results are consistent over the ridge 
itself. 

The study demonstrates the importance of the lower atmosphere in mountain 
dynamics. The vertical mixing of the PBL on sloping terrain produces unique effects and 
it is clear that the inclusion of a realistic PBL parameterization is even more critical in 
these frontogenesis simulations than in their flat topography counterparts. To complete 
the description of the frontogenetical forcing, the terms represented by equation (45) must 
be exphcidy calculated. The vertical mixing effect depends heavily on the initial potential 
temperature distribution and future studies should also address this sensitivity. In 
particular, the inclusion of diurnal surface heating should have a major effect of frontal 
structure and dynamics. Also, the initial temperature state is not realistic in the PBL 
simulations, and the adjustment while the PBL forms may introduce extraneous frontal 
effects in the first stages of the integration. Similarly, the sensitivity to the boundary 
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layer constants, particularly the mixing length / should be addressed. Ultimately, the 
experiment should be expanded to a three-dimensional domain to evaluate the PBL effects 
on the along-ridge variations found by Gross (1994). 
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